Spin-spin correlations of entangled qubit pairs in 
the Bohm interpretation of quantum mechanics 



A. Ramsak 

Faculty of Mathematics and Physics, University of Ljubljana, Ljubljana, Slovenia 
J. Stefan Institute, Ljubljana, Slovenia 

E-mail: anton . ramsakOf mf . uni-1 j . si 

Abstract. A general entangled qubit pair is analyzed in the de Broglie-Bohm 
formalism corresponding to two spin-1/2 quantum rotors. Several spin-spin correlators 
of Bohm's hidden variables are analyzed numerically and a detailed comparison with 
results obtained by standard quantum mechanics is outlined. In addition to various 
expectation values the Bohm interpretation allows also a study of the corresponding 
probability distributions, which enables a novel understanding of entangled qubit 
dynamics. In particular, it is shown how the angular momenta of two qubits in this 
formalism can be viewed geometrically and characterized by their relative angles. For 
perfectly entangled pairs, for example, a compelling picture is given, where the qubits 
exhibit a unison precession making a constant angle between their angular momenta. 
It is also demonstrated that the properties of standard quantum mechanical spin- 
spin correlators responsible for the violation of Bell's inequalities are identical to their 
counterparts emerging from the probability distributions obtained by the Bohmian 
approach. 
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1. Introduction 

It is well known that several phenomena pertinent to spin-| systems can not be described 
in terms of classical variables although to some extent a geometric terminology may be 
applied to a spin state of a single electron which is customary presented on the Bloch 
sphere, where the polar and azimuthal angles can be viewed as the Euler angles of a 
unit vector pointing along the spin direction. Still, within the framework of standard 
quantum mechanics there is no answer to the questions like: How could one imagine spin 
properties of an entangled electron pair in terms of classical variables? The question is 
related to the problem of quantum entanglement, which is a rather mysterious concept 
and has remained so since it was first introduced by Schrodinger back in 1935. Despite 
this, with the recent surge in interest stimulated by quantum information processing 
[HE] an d quantum sensing below the quantum limit [3], discussions about the meaning 
of entanglement and how it might be visualized beyond the standard quantum mechanics 
have once more come to the fore [H H]. In this context can also be set recent intense 
exchange of views whether the wavefunction is physically real and can be directly 
measured or represents only a statistical tool that reflects our ignorance of the particles 
being measured [6]. 

Quantum entanglement is strongly related to the discussion of locality in quantum 
mechanics and with the introduction of Bell's inequality also to the question of the 
existence of hidden variables [7J [8] . In order to set experimentally observable boundaries 
to various hidden variable theories several other inequalities have been introduced 
[9j [10]. The outcome of numerous experimental tests of such inequalities in the last 
three decades performed on two-particle systems ruled out all local hidden variable 
approaches, although it seems that non-locality concept itself is not sufficient to be 
consistent with quantum experiments [UJ . However, for multipartite entangled states 
local hidden variable models can be introduced such that the violation of most of the 
Bell- type inequalities does not imply and new inequalities have to be established [HI [13] . 

In view of the fact that quantum entanglement is not an observable according to 
the usual rules of quantum mechanics, it does not have a direct classical analogue. 
However, some insight can be given by a description based on the pilot-wave approach 
of de Broglie [2] and later developed by Bohm [15] - where an effective non-local 
potential is introduced along with hidden variables. This theory is manifestly non- 
local, i.e., non-local potentials represent non-local action which does not get small with 
spatial distance between the particles and this is the key ingredient which guarantees 
predictions in accordance with standard quantum mechanics. While it seems that such 
an introduction of non-local variables does not lead to new predictions - it might turn 
out to be interesting in the context of possible device independent analysis of entangled 
systems [161 HB]- In particular, one can ask questions related to pre-existing properties 
of the system; explicitly, what each individual particle in the observed system is doing 
prior performing the measurement. 

Since the introduction of the Bohm interpretation of quantum mechanics it has 
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been in detail explored how this approach leads to outcome predictions of various 
experiments identical to the results of standard quantum mechanics, and gives due 
to its nature, additionally also the interpretation of the motion of particles in terms 
of particular coordinates - commonly termed as the Bohm hidden variables. After the 
establishment, the Bohm formalism was extended from original single spinless particles 
to many-body fermionic or bosonic systems pTTj [TBI IT9] and also to the quantum field 
theory, including creation and annihilation of particles [201 121] ■ Recently it was shown 
how the Born-rule probability densities of non-relativistic quantum mechanics emerge 
naturally from the particle dynamics of de Broglie-Bohm theory [22] • The interpretation 
of several typically quantum experiments was developed, for example, for the double slit 
experiments or tunneling of particles through barriers [T7], the Aharonov-Bohm effect 
[23J, and Stern-Gerlach experiments [24]. The Bohmian theory in terms of well-defined 
individual particle trajectories with continuously variable spin vectors was successfully 
applied also to the problem of Einstein- Podolsky- Rosen spin correlations [25]. Among 
the simplest quantum objects are two state systems and the formal Bohm approach 
is consistently extended to spin-| systems with non-relativistic formalism based on 
the Pauli equation [25], causal rigid rotor theory [27], via spinor wave functions [18J, 
the Bohm-Dirac model for entangled electrons [28J, by Clifford algebra approach to 
Schrodinger [29J or relativistic Dirac particles [30J. 

In this paper we concentrate on the question of the visualization of the angular 
momenta - spins - of a system of two entangled particles by the approach of the causal 
interpretation of a system of two spin-| particles as introduced by Holland [27J,[T7] where 
the starting point is the mapping between the quantum rigid rotor and a spinning top in 
the presence of a quantum potential. In particular, we apply this ontological formalism 
to the case of a two-particle state 

l*> = cos t |n) + e* sin t ||t> , (1) 

where for the sake of convenience we use spin-| notation where |n) corresponds to 
the state of the system when the first particle (qubit) is in the "up" state, i.e., in the 
direction of the z-axis and the second qubit is in the state "down". The generalization 
to states spanned by |tt) an d or to systems represented as mixed states is possible, 
but it is not considered here. Qubits are not restricted to real spin of electrons and may 
be realized by any appropriate two state quantum system, as are for example flux qubits 
in superconducting rings [31], charge pseudo-spin of electron pairs in double quantum 
dots [32], two flying qubits in quantum point contacts [33], entangled photon pairs [34], 
or even most recently studied two-qubit composite systems ]35j. 

The causal approach considered here has been introduced more than two decades 
ago but exact solutions are known only for the case of a single qubit and for some 
properties of a qubit pair in the singlet state. It is also known that the Bell inequalities 
are violated in this formalism and recently the quantum entanglement of a qubit pair 
was linked to a particular motion of angular momenta in the Bohmian space of hidden 
variables in a way which is entirely missed in the usual quantum mechanical approach 
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[3T] . The problem of the motion of two entangled rotors in the Bohmian space of 
hidden variables can be efficiently treated numerically, but the author is not aware of any 
quantitative investigation of this problem. The aim of the present paper is to fill the gap 
by a comprehensive review of spin-spin correlations and the corresponding probability 
distributions of Bohmian hidden variables in a system of two entangled qubits in the 
quantum state Eq. ([I]). We explore also the relation to the Bell's inequalities for the 
case of a perfectly entangled qubit pair. 

The paper is organized as follows: after the introduction in Section 2 the model for 
two quantum rigid rotors in the causal picture is introduced. In Section 3 numerically 
extracted probability distributions of various hidden variables related to two spin-| 
rotors and the corresponding ensemble average values are presented. The Shannon 
entropy for a discretized probability distribution of the angular momentum 2-axis 
projection is given. Section 4 is devoted to geometrical aspects of two quantum rotors, 
in particular, to the study of the relative angles between the angular momenta and the 
corresponding fluctuations. In Section 5 we express the Bell's inequalities in terms of 
the correlators derived in Section 3. Section 6 is devoted to summary and conclusions. 



2. Bohmian treatment of rigid rotors 

2.1. A single rotor 

As a preliminary to the detailed examination of entangled rotor pairs we consider a 
single quantum spherically symmetric rigid rotor treated by the Hamiltonian 

M 2 

H = w , (2) 

where M represents the angular momentum operator - a differential operator whose 
components are the infinitesimal generators of the rotation group SO (3) and / is the 




Figure 1. (a) Rotor 1 principal axis ei with Euler angles Ai = {o:i,/3i,7i} (e2 for 
rotor 2 is not shown) and angular momenta Mi^, for a general case of two qubits. In 
the case of unentangled rotors, z? = 0, the vectors Mi, ei and the z-axis are coplanar, 
Eq. (|19[) . Euler angles 9%, <f>i correspond to the angular momentum Mi and #2, 4>2 (not 
shown) to M2. (b) Two qubit case with shown relative angle $ and the xy- in-plane 
azimuthal angle <j> between Mi and M2. 
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moment of inertia [T7J |27] . Eigenstates ip of the three mutually commuting operators 
M 2 , M z and e ■ M are functions of Euler angles A = {a, (3, 7}, specifying the orientation 
of a rigid body with the principal axis defined by a normalized vector e. The angular 
momentum eigenvalues are defined by 

MV =s(s + l)ip, (3) 
M z ip = mil), (4) 
e ■ Mip = nip, (5) 

where s is the total angular momentum quantum number, and m, n are the angular 
momentum projection quantum numbers: in the spatial z-axis direction and with respect 
to the rotor axis e, respectively. 

In order to analyze ontological properties of rigid rotors in de Broglie-Bohm space of 
hidden variables, we in this Section closely follow the formalism introduced by Holland 
[T71 127] . where a quantum rigid SU{2) rotor is represented by an ensemble of SO (3) 
rotors. The standard quantum mechanics does not describe coordinates and trajectories 
of particles, but only expectation values of various observables and the probabilities for 
particular outcomes of experiments. On the other hand, in the Bohmian interpretation 
of quantum mechanics, one can additionally discuss positions and momenta of particles 
and the probability densities that particles actually are at particular positions with a 
particular momentum, not only the probabilities that such positions or momenta will 
be obtained after completing the measurement process. 

Following Bohm the wave function is expressed as 

4>(a,P, 1 ) = Ke iS , (6) 

where 1Z and S are real functions [15J. In this approach the angular momentum in the 
Bohmian space is given by a real, three dimensional vector M = iM5, 

M x = - cos f3 dS/da + sin cot a dS/d(3 - dS/d^ } (7) 

sin a 

M v = sin p dS/da + cos f3 cot a dS/d/3 - dS/d-f, (8) 

sin a 

M z = - OS I '3(3. (9) 

The dynamics is determined from the Hamilton- Jacobi-type equations for the classical 
Hamiltonian with an additional quantum potential Q, 

H=— + Q 1 Q = -^. (10) 

The quantum potential Q generates the quantum torque T = — WVQ acting on the 
rotor as dNL[\(t)]/dt = T[A(t)]. The equations of motion reduce to a set of first order 
non-linear differential equations for the trajectories in the configuration space, 

la = dS/da, (11) 
1$ = (dS/d(3 - cos a dS/dy)/ sin 2 a, (12) 
J 7 = (OS/d^ - cos a dS/d/3)/ sin 2 a. (13) 
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where particular trajectories X(t) in the space of hidden variables are the ensemble 
representatives, determined by the initial condition A(0) = {a , (3 , 7 }. 

2.2. Spm-\ 

In the analysis of spin s = | systems we chose solutions with spin projection along the 
rotor axis n — ^ an d along the z-axis m = ±|, where the corresponding eigenstates are 
the Wigner matrices D^, 

Uf(X) = — -!=!?! (B, a, 7) = ^Le" i/3/2 - in7 cos -, (14) 

= -^=n\ n (P, «> 7) = sin f • ( 15 ) 

One should note that the core of this approach lies in the coexistence of the motion of a 
rigid top whose configuration space is 5*0(3) but its motion is guided by a wave and the 
quantum potential whose spin configuration space is SU(2), as discussed in Ref. jTTj- 
Wave functions ip are normalized, and 1Z 2 represents the probability density for A, with 
the normalization condition in the domain A 

III n 2 smadad/3d-f= TZ 2 dX = l. (16) 

JO Jo Jo J A 

The quantum equilibrium ensemble average of some function fi(X) is then given by 

(//) = / fi(X)n 2 (X)dX. (17) 



A 



As an example, consider a qubit in the "spin up" state, tj) — (A \\) = u^. The exact 
solution is given by X(t) = {a , f3 — t/T,j — t/r}, where {«o,/3o;7o} are the initial 
values for A and r = 41 cos 2 (ao/2). The angular momentum vector precesses uniformly 
anticlockwise about the z-axis with a constant polar angle 9 = ao/2, 

M(t) = |M| (sin 9 sin 0,sin0 cos 0, cos 9), (18) 

|M| =tt\, <P = Po-t/r. (19) 
2 cos 9 

Note that M z and e ■ M are constant |, while the angular momentum vector length is 
constant, |M| > |, dependent of the initial condition for X{t). 

2. 3. Entangled qubit pair 

The generalization to the case of two rotors is straightforward [17J. We consider here 
a general two qubit state with vanishing total angular momentum projection, Eq. ([TJ). 
The guiding wave function ip(X) = (A|^), 

V(A) = cos ^(AiHCAa) + sin ^(X^u^) , (20) 

is given in six dimensional space spanned by A = {Ax,A 2 }, with Ai 2 representing the 
coordinates of the first and the second rotor, respectively. The angular momenta are 
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determined by the set of equations (jMH]) and (I11H13I) . generalized to the case of two 
qubits, M 12 = iM 1)2 iS. The corresponding Hamiltonian is given by 

g _ Mf + Mj (Ml + Ml)K 

21 2IK ' 1 ' 

It should be noted, that even for two non-interacting, but entangled qubits, the quantum 
potential Q(X) - the second term in Eq. ( 12T|) - represents an instant interaction between 
the rotors as a fingerprint of the quantum nature of the problem. The solutions for each 
of the angular momentum vectors Mi and M 2 are functions of six common coordinates 
forming the trajectory X(t) determined by six initial values A(0). The ensemble average 
is for the case of two qubits also given by Eq. f lTT]) . but in six-dimensional space and 
with dA = dAidA2- 

For the case of non-entangled qubits, $ = 0, the solutions are given by appropriately 
applied Eq. ( Fl9l) . while for a general $ the angular momenta and the rotor axes are not 
coplanar, Fig. [TJ(a), and exhibit a rich variety of precessional motions. Total energy 
E = s(s + 1)// is a constant of motion equal for all A. Even though also some other 
constants of motion can readily be identified and one can prove that the orbits X(t) are 
periodic in particular subspaces of the phase space, the time evolution along the orbits 
can not be tracked analytically in general. The analysis of the topology of the orbits 
reveals several classes forming distinct portions of the phase space and the survey of 
these properties of the model will be presented elsewhere [36J. 

3. The probability distributions 

Here we focus on the analysis of various static, t = 0, probability distributions of 
two spin-| rotors and we compare particular ensemble averages with the corresponding 
results obtained by standard quantum mechanics. All probability distributions and the 
ensemble averages considered here are time-independent while some quantities specific 
for time- dependent properties of the rotors are presented in Ref. [37] . 

3.1. Single particle properties 

First we consider single particle properties of the rotor pair. For the state with 
a vanishing total spin projection the equality Mi z {\) + M 22 (A) = is fulfilled for all 
A and the projections along the rotor axes are constants of motion by construction, 
2 ■ Mi ;2 = |. The average values of the angular momentum calculated by Eq. (fTTj) 
coincide with the standard result for spin expectation values 

(Mi) = -(M 2 ) = ^(0,0,costf) = (tf|Si|tf>, (22) 

where Si = \(o~i X) o~i y , o"i 2 ) is the ordinary spin operator for the first particle and a are 
the Pauli matrices. 

Despite that average Bohmian momentum Mi and standard quantum counterpart 
Si are equal in this case, there is a significant difference of both angular momentum 
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is 



vectors: the state |\&) is an eigenstate of S^, which is dispersionless, while |Mi 
different for different members of the ensemble and distributed according to some 
probability distribution corresponding to the relation Eq. ( fl~9l) . In order to gain the 
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Figure 2. The probability distribution dPjd[i of /i = Mi^ for various i? (blue lines). 
For unentangled qubits it is given by dP/d/j, = 1(1 + 4/x 2 )~ 5 / 2 while the distribution 
for /1 = Mi z (red lines) is for $ = a delta- function at jit = — |. The distributions are 
independent of if and are shifted vertically for clarity. For /1 = |Mi| 2 the distribution 
is given analytically by dPjd[i = 0(/i — j)(2fi)~ 3 [derived from Eq. (|26l) ] and is 
independent of 1? (black line). The probability distribution of the rry-plane projection 

of Mi, n = 



M xy , is shown with green dashed lines. Thick vertical 



lines represent (M2 Z ) (left set) and (M xy ) (right set). 



insight into the properties of the Bohmian motion of the rotors we analyse this and 
various other probability distributions corresponding to hidden variables. In general 
the quantum equilibrium probability distributions corresponding to various functions 
//(A) studied here are defined by 

d ^tL = f S[^-^\)}n 2 (\)d\, (23) 

where 5(fi) is the Dirac delta function. Average values corresponding to such 

distributions are then expressed, e.g., for //(A) = M lx , by 

r r°° dp 

(M lx ) = / M lx (A)ft 2 (A)dA = / M lx — — alM lx . (24) 
Ja J-oc dM lx 
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Therefore dP represents the probability that for a given qubit pair in the state |\?) 
the x-axis angular momentum component of the first particle is within the interval 
[M\ x , M\ x + dMi x ] not the probability that Mi x will be obtained after the performed 
measurement. Note that the two probability densities are in general not equal, while 
the probability for the outcome of a particular experiment is always identical for both 
approaches, the Bohm and the standard quantum mechanics. 

For few limiting cases the probability distributions can be derived analytically, while 
in general one has to rely on a numerical evaluation of Eq. ( 1231) . for example, by the 
limit of a sum over the configuration space, 

dP< ) i=N i=N 

"^T = fe E^-MA*)]^)/^ 2 ^), (25) 



+ 



where 8 e (/i) is a normalized rectangular function of width e and iV is the total 
number of samples Aj forming a uniform grid in the configuration space spanned 
by {cos oti, 0i, 71, cos ^2, fa, 72}- The hidden variables corresponding to the internal 
rotation around the principal axes of the rotors are irrelevant for the quantities 
considered here, thus we fixed 7^2 = 0. To obtain sufficiently converged results with 
Eq. (125]) we used N = 1024 4 and e = 10~ 3 for the full domain A corresponding to 
I cos ai^ I < 1 and < fa$ < 2n while for particular cases with higher symmetry the 
domain can be reduced. The error due to the finite grid used in the calculations presented 
throughout the paper is always below the width of the lines in the corresponding figures. 

The probability distribution of momentum length dP/d\Nli\ and the corresponding 
moments can be for i? = derived analytically, 



dp _ e(|Mi| - i; 



d|Mi| 4|Mi| 



(26) 



(IM^ 2 ) = i = !<¥|S?|¥), (27) 

where 0(/x) is the Heaviside step function. 

All single particle distributions are independent of (p, and the numerical analysis 
gives a firm evidence that the probability distribution Eq. fl26|) is in fact also independent 
of Numerically we determined also the probability distribution for the xy-plane 
projection of the angular momentum, /i = yj M\ x + M\ y = M xy , for unentangled qubits 
given analytically by dPjd^i = 16/x(l + 4/x 2 ) -3 . The distribution is shown in Fig. [2] for 
various $ and the average (M xy ), right set of thick vertical lines in Fig. |2l ranges from 
f to § from the unentangled to the fully entangled case, respectively, which is less than 



the standard quantum mechanics counterpart (\&| yj Sf x + Sf y \^/) = 

Additionally are shown the probability distributions for \x = M\ x and fi = M2 Z = 
—Mi z . The average (M2 Z ) = —^cosfi is shown by the left set of thick vertical lines. 
The corresponding probability distribution progressively develops from a delta-function 
for unentangled qubits into a distinctive probability distribution curve, constant for 
I A* I < I' ^ or entanglement at d = tt/2. The characteristic kinks in the distributions 
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at = | have the origin in Eq. §5§ with the eigenvalue n = ^ which sets the lower 
limit for the magnitude of the angular momentum in Eq. fl26|) . 

For perfectly entangled qubit pairs the probability distribution of the polar angle 
cos#i = M lz /|Mi| [Fig. []Ja)] is constant, dP/dcos9i = ~, as anticipated for a singlet 
(but valid for any <p, as argued above). As a direct consequence emerge the exact 
expressions for the probability distributions of \i = M\ z , \i = M\ z (identical are the 
results corresponding to \i = M lx ( y )) and fi = M xy , 
dP 4 

— = -min(l,|2 /U |- 5 ), for^ = M lz , (28) 
dfj, 5 

dP 8 

— = -min[(4/i)- 1 / 2 ,(4 / i)- 3 ], for fj, = M\ z , (29) 

^ = T^t 1 - 6 ^ - ^ + 2 ^ + 6/)^- V], for /i = M xy . (30) 

Finally, let us comment on an interesting property following from Eq. f|27j) . The 
average square of the Bohmian momentum represents exactly | of the corresponding 
quantum value, which can be by taking into account the equality of the Bohmian average 
energy and the quantum expectation value, E = (H) = (^\H\^/), re-expressed into the 
identity 

'm 2 12 h\ 1 



n / 2 



<|Mi, 2 | 2 >. (31) 



The identity represents the virial theorem: for an arbitrarily entangled, but non- 
interacting, two qubit state the ensemble average of the quantum potential in Eq. fl2Tl) 
represents exactly one half of the average kinetic energy in the Bohm space of hidden 
variables, 

(Q) = i(H-Q). (32) 

This means that the quantum potential contribution is never negligible, (Q) = E/3, 
which indicates the contextual nature of considered hidden-variables formulation of 
quantum mechanics [8], i.e., various angular momentum correlation functions after the 
performed measurement are in general different from Bohmian pre-measurement values 
where the rotational energy is diminished on the account of the potential part. This is 
similar to the more familiar case of the orbital motion of point-like particles where the 
Bohmian particle velocity correlations differ from the canonical momentum quantum 
counterparts [38] . 

3.2. Entanglement of formation and the Shannon entropy 

Apart from the ontological point of view, an important distinction between standard 
and de Broglie-Bohm interpretation of quantum mechanics reflects in the probability 
distributions which can be discrete in one case and continuous in the other [39]. Due to 
this qualitatively different nature of the probabilities the entanglement of a qubit pair 
in the two approaches is quantified differently 
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The probability distribution for a measurement outcome of spin projection Si z (or 
M\z) is discrete with the probabilities pi = cos 2 | and P2 = 1 — pi for "spin up" and 
"spin down", respectively. An appropriate measure for the degree of entanglement is the 
entanglement of formation [JO], the asymptotic conversion rate to maximally entangled 
states from an ensemble of copies of a non-maximally entangled state [2]. It is given by 
the von Neumann entropy S(pi) = — Trp x log pi, where pi is the reduced density matrix 
obtained by tracing the qubit pair density matrix over the degrees of freedom of particle 
2. The log is to the base 2 and entropy is expressed in bits. In our case of a pure state 
with p = the von Neumann entropy is identical to the binary (Shannon) entropy 

H s = ^ log ^, (33) 



where is i = 1, 2 for S\ z — ±~, respectively. 

In de Broglie-Bohm approach p = M\ z is a continuous random variable thus the 
probability distribution p(p) = dP/dp is not discrete, as discussed in the previous 
subsection. In this case an analogue measure of entanglement is the Shannon entropy 
H± for a random binary variable that the angular momentum vector of the first particle 
is in the upper or lower hemisphere with probabilities p + = J °° p(p)dp and p_ = 1 — p + , 
respectively. Although (M\ z ) = (Si z ), Eq. ( 1221) . the probabilities pi and p + are in 
general not equal as is shown in Fig. [3] (dotted line), where are as a function of pi shown 
also the binary entropy H± (dashed) and the entanglement of formation (thick full line). 




cos 



Figure 3. Probability p+ for Mi z > as a function of pi — cos ^ (dotted line, 
right scale). Note the deviation from a straight line which reflects in the inequality of 
the binary entropy H± (dashed) and the entanglement of formation S(pi) (thick full 
line). Thin lines represent the renormalized Shannon entropy H u /v corresponding to 
bin widths A = 2~ v = \, i 



Since the probability distribution p(M l2 ) is a continuous function an appropriate 
measure is given also by the differential entropy, 

h = - / p(p) logp(p)dp. (34) 
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The differential entropy is related to the Shannon entropy Eq. (1551) with quantized 
probability = p(/ij)A, where the range of [i is divided into bins of length A and \ii 
is a value of /i within the i-th bin such that pi = f^i^ p(n)d[i [H]. For the sake 
of convenience we take A = 2~ v which leads to a relation valid for the corresponding 
Shannon entropy H v , 

H u — v h, as v — t- oo. (35) 

This means that v + h is approximately the number of bits on the average required to 
describe M\ z to v bit accuracy. The entropy H v is zero for unentangled qubits and is 
maximum for a fully entangled case. Maximum of h, calculated from Eq. (128]) . equals 
hmax — log f e 1 / 4 ~ 0.68 with the corresponding support length containing most of the 
probability Vol(Mi2) = 2 hmax pa 1.61, which can also be estimated visually from Fig. |2j 
In Fig. |3jis with thin lines presented normalized Shannon entropy H v fv for v — 2, 4, 8. 
We found an excellent agreement with the asymptotic relation Eq. (135]) for v > 6 where 
the deviations are below the width of the lines. 

The probability distribution dPjdM\ z thus determines the degree of entanglement 
similarly to the conventional approach if binary entropy criterion is used. Quantitative 
differences between H± and S(pi) vanish for unentangled or fully entangled pairs, as 
expected. The main difference between de Broglie-Bohm and the standard view is, 
however, that the information contained in de Broglie-Bohm probability distribution is 
due to its continuous nature dependent on the discretization, therefore can be arbitrarily 
large, ~ v + h. It should be noted that all this information is not accessible and 
probably can not be used for a transfer of information because during the process of 
spin measurement the spin projection outcomes will become ±| - i.e., a non-continuous 
value - thus the entanglement of formation is limited to at most one bit per entangled 
qubit pair. 



3.3. Two-particle distributions 

One of the main motivations of the present paper is the analysis of the two-particle 
properties which reflect quantum mechanical correlations of two entangled qubits. We 
consider here the probability distributions dP/dfi for /i = M lx M 2x and /i = M lz M 2z . 
The probability distributions are presented in Fig. H]for various degrees of entanglement 
and for (p = 0. 

The distribution of M\ z M 2z is independent of (p and is for two unentangled qubits a 
delta function at jj, = — | while for fully entangled triplet state it takes a form identical 
to the distribution of M lx M 2x but with opposite sign of fi. For the singlet the two 
distributions are identical (not shown). Due to the equality M\ z = —M 2z one can for 
d = n/2 from Eq. (129]) express these distributions analytically, 
dP 8 

— = -min[(4r / /i)- 1 / 2 ,(4r / /i)- 3 ], for /x = M lz M 2z , (36) 
djj, 5 

dP 

— = (47//1)- 1 / 2 , for fj, = M lz M 2 ,/(|M 1 ||M 2 |), (37) 
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Figure 4. Probability distributions dP/dfi for /j, = M\ x Mi x (blue lines) and 
fi = M\ z Mi z (red lines) for various i? and y> = 0. The distributions are vertically 
shifted for clarity. Thick vertical lines represent (Mi z M2 Z ) and (M\ x Mi x ) (left and 
right set, respectively). 



where 77 = — 1, while the corresponding distributions of the x- and ^/-components of 
the angular momenta have the same form for the singlet or the triplet, with 77 = =pl, 
respectively. Here the kinks in the distributions are at = \. Note that /i _1//2 in 
Eq. (1371) reflects the isotropic distribution of the angular momentum components. Let 
us also point out that for every A the angular momentum vectors have equal length, 
I Mi (A) I = |M 2 (A)|, while for ■& < f this equality holds only on average, (|Mi|) = (|M 2 |), 

I Mi I — |M 2 | is finite, for example, 



i.e., the width of the probability distribution for fi 
(/i 2 ) = I is maximum at 1? = 0. 

The ensemble average values (Mi z M2 Z ) and (Mi x M2 X ) are indicated in Fig. H] by 
thick vertical lines. In general, the ^-dependence of correlators (MuM 2 j) is identical to 
the standard quantum mechanical spin-spin correlation tensor, 

' sin $ cos ip — sin $ sin <p ' 
sin $ sin (p sin t? cos <p 
-1 



(*\S U S. 



2/ 



1 4 



V 



(3* 



/ 



where represent the corresponding x, y, z components. In Fig.^a) are shown various 
spin-spin correlators {MuM 2 j) as a function of 1?. The dependence of i? is slightly 
different in comparison with Eq. (j38p but at 1? = | the Bohmian correlators are exactly 
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| of their quantum counterparts, as is the case for the single-particle averages discussed 
in the previous subsection, Eq. ( |2"7|) . The second moments (Mf z ) and (Mf x ) are also 
shown for a comparison with standard quantum mechanics result (Sf xy z ) = \. 




Figure 5. (a) Renormalized spin-spin correlators: (M\ X M2 X ) I cos ip or 
(M ly M 2x ) / sin ip (thick full line). (Mi x M 2x ) / cosfl cos ip (dotted line), and <p- 
independent (M lz M 2z ) = ~{M? Z ) (dashed-dotted line) and {M? x ) = (Mf y ) (dashed 
line), (b) Angular momentum correlator (Mi • M2) (full lines) in comparison with 
standard quantum mechanics result (Si • S2) (dashed lines) for two values ip = and 
ip = 7T. 



In Fig. E^b) is with full lines presented the correlator (Mi ■ M 2 ) for singlet and 
triplet-type of qubit pair configuration in comparison with quantum mechanical value 
given by (^\S\ ■ S2 1 1 ^ r ) = |(2sin , (9cos(/? — 1) [42J. For unentangled state all correlators 
are equal to —4 while for = 7r/2 the Bohmian correlators for any f exhibit anticipated 
relation (Mi ■ M 2 ) = |(*|Si • S 2 |*) which together with Eq. (ETJ) leads to the identity 

((Mi + M 2 ) 2 ) = |<*|(S X + S 2 ) 2 |vI/>. (39) 
4. Geometrical view of a qubit pair 

In standard quantum mechanics a visualization of a relative motion of two angular 
momenta of an entangled qubit pair is not given directly, although by appropriately 
defined angle operators some imagery can be given [12]. On the other hand, in de 
Broglie-Bohm interpretation of quantum mechanics an analysis of geometrical quantities 
as are, for example, relative angles between the momenta in the space of hidden variables 
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is simple and directly feasible. We apply the approach introduced in the previous 
Section and systematically present geometric quantities relevant to the qubit pair given 
by Eq. (Op). We consider two possibilities: the angle made by the angular momenta and 
the azimuthal angle made by the x?/-plane projections of the momenta, as follows. 



4-1. Angle between two angular momenta 

The angular momenta Mi and M 2 of two qubits make an angle $ in the Bohmian space 
of hidden variables A as shown in Fig. [T](b). The average cos$ and the corresponding 
dispersion A cos $ are given by the expressions 

(cos4) =/ff,ff,\, (40) 
x 1 \|Mi||M 2 |/ 

(A cos $) 2 = (cos 2 $) - (cos $) 2 . (41) 

The average cosine is presented in Fig. |6](a) as a function of °d and for various ip. 
For i? = 7r/2 we find exact relation (cos$) = |(2cosy? — 1), identical to the standard 
quantum mechanics expression [12]. As expected, perfect anti-parallel alignment is 
found for the singlet state {ip = ir, red line), while momenta for the triplet state, ip = 0, 
are only partially aligned, (cos$) = |. In Fig. [6](b) is shown the variance Acos$. 
For (p ~ 7T it is strongly suppressed for a wide range of $ > | indicating a significant 
stability of anti-parallel angular momenta configuration of the singlet qubit pairs (red 
line). The variance for states close to the triplet configuration, on the other hand, is 
much higher and is of the order of average cos $. 




Figure 6. (a) (cos $) vs. $ for various (p. Bullets represent exactly determined special 
cases, (b) Variances Acos$ corresponding to (a). 



The correlators applied in Eq. PU|4T]) exhibit a similar symmetry structure as 
(y\S u S 2j \V) in Eq. fl38D, 



M u M 2j 



1 

3 



B x cos ip —B x sm.(p 
B x sin ip B x cos ip 
-B, 



\ 



(42) 
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where B xz are functions of shown in Fig. [7Ja,b) (red lines) and (cos$) = 
\[2B x {$)cosip-B z {-&)\. 

For comparison is in Fig. UJa,b) additionally shown the Bohmian counterpart of 
sin$ in Eq. (ESJ), C M = 6(M lx M 2x ) / cosy? = 6 (M ly M 2x ) / sin <p (dotted line). Let us 
remark that both, B x and Cm, are very close, though not identical to sin-!?. 




Figure 7. sini? (thin line), Cm (dotted line), Cb (equal to (cos0) for ip = 0, dashed 
line), and B x {d) (red thick full line). Variances (Acos</>) 2 and (Asint/)) 2 are shown 
for ip = or 7r (dashed-dotted lines). For ip = ±ir/2 the two curves are interchanged. 
B z (i9) (red long dashed line) and the same set of quantities as in (a) , but renormalized 
by sin i?. Note the identity of all curves at $ = ^. 



4-2. Relative azimuthal angle 

We consider here qubit pairs with vanishing total z-axis projection of spin, therefore 
the angles between projections of My onto the xy-plane play a special role [37]. We 
analyze the difference of azimuthal angles of the angular momenta <f) = <p2 — 4>i via 
properties of average cosine and the corresponding variance defined by 

, , x / M lx M 2x + M ly M 2y 

(cos0) = ( = 

\^(Ml + Mf y )(Ml + Ml y ) 

(A cos (pf = (cos 2 <p) - (cos <p) 2 . (44) 

By inspection of the wave function ip, Eq. ( 1201) . we find that a finite ip represents only 
a shift of one of the azimuthal angles for each member of the ensemble, /?2 — > /?2 + ip, 
which results in the identity (0) = ip and the decoupling (cos0) = (7b cosy?, where 
Cb is a function of •& only. An analogous result holds for (sin <fi) = Cb sin ip, thus 
(cos(0 — ip)) = Cb and (sin(0 — ip)) = 0. In Fig. [TJa) is Cb($) presented by a dashed 
line in comparison with sin$ (thin full line). The azimuthal angle fluctuations can 
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be analyzed by variances A cos and A sin 0, also shown in Fig. [7J(a) (dashed-dotted 
lines) . Note that with increasing d suppressed angle fluctuations signal a higher degree of 
entanglement as a highly correlated distribution of angular momenta making azimuthal 
angles progressively closer to (p. 



5. Bell's inequalities for a Bohmian ensemble 

One of the most striking properties of two entangled qubits is the violation of Bell's 
inequality [7] derived for systems described by local deterministic hidden-variable 
theories. For more general versions of the inequality, Clauser, Home, Shimony, and 
Holt [13] showed that 

|P(a, b) + P(a, b') + P(a', b) - P(a', b') | < 2, (45) 

where P is the counting correlation of measured outcomes of spin projections along 
polarizer directions given by normalized vectors a, b, a' and b' [7J. It is readily seen 
that for the quantum prediction 

P(a, b) = (singlet|(a • c?i)(b • <7 2 )|singlet) = —a • b (46) 

with well chosen polarizer directions the inequality can be violated by a factor as large 
as \/2. 

Theory of two entangled qubits presented here is based on the hidden-variable space 
spanned by A, but the formalism is non-local due to the quantum potential generating 
instantaneous interaction between the rotors. The question arises then to which degree 
the inequality Eq. (T4"5l) is respected for the case of two rotors represented as a Bohmian 
ensemble. 

Let us first reconsider spin-spin correlations C(a, b) = (\l/|(a • e?i)(b • cr 2 ) | ^ ) for a 
general state Eq. ([I]). Correlations can be for convenience expressed directly in 
terms of spin operators Si i2 and Sf 2 as, 

C(a,b) = 3(*|^^)|*)= (47) 
= (o>xb x + a yb y ) sin i? cos <p + (a y b x — a x b y ) sin i? sin <p — a z b z . 



In a similar manner one can construct for the case of Bohmian rotors an analogous 
correlator, related to the quantities B x z studied in the previous section, 

»b).3/' a 'S, M, 'U (48) 



\ |Mi[|M 2 | / 
a x b x + a y b y )B x cos y? + (a y b x - a x b y )B x sin <p - a z b z B z 



For a perfectly entangled qubit pair 



|¥) = -L(|t|) + e*|;t», (49) 
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the correlations C(a, b) and B(&, b) are identical due to the equality B x = B z = 1, see 
Fig. UKjo) and also Eq. (|37j). 

Bohmian correlation measure defined by Pb = B(&, b) thus for any setup of 
polarizators violates the Bell inequality Eq. ( H5|) exactly as the quantum counterpart 
P does. It is well known that Bohmian and standard approach of quantum mechanics 
are isomorphic regarding the prediction of experimental outcomes therefore the equality 
of spin-spin correlations in both approaches apparently is not a surprise. However, 
one should remark, that the inequality Eq. ( H5l) is derived under the assumption of 
discreteness of the outcome of polarization measurements, while in the present Bohmian 
analysis the spin of rotors and the corresponding projections can take any value. More 
importantly, Pb corresponds to instant correlations of an undisturbed rotor pair, prior to 
accomplished measurement, not to the probabilities for particular outcomes of discrete 
values obtained after the measurement of spin (such probabilities are trivially identical 
in all approaches). Therefore, although the correlators of the continuously distributed 
Bohmian angular momenta are identical to the quantum counterparts, those momenta 
are not directly related to the results of quantum measurements of spin projections of 
individual qubits, but during the process of the polarization measurement the entangled 
rotors adiabatically develop the measured discrete spin values from pre-measurement 
distributions of continuous angular momenta. This result resembles the findings of 
an alternative analysis of Stern-Gerlach experiments or Einstein-Podolsky-Rosen spin 
correlations where de Broglie-Bohm causal predictions of experimental outcomes are 
identical to standard quantum mechanics results (2H [25] . 

6. Summary 

We concentrated on time independent properties of a two spin-| particle system, as 
viewed in the Bohmian space of hidden variables. The state Eq. ([1]) represents one 
of the simplest many-particle quantum systems exhibiting quantum entanglement and 
although it is determined by two parameters only, the angular momenta in the Bohmian 
space display an infinite range of detail which we endeavoured to present by means of 
various probability distributions and the corresponding average values for the Bohmian 
ensemble. The main motivation was the analysis of entanglement aspects of a qubit 
pair, in particular the exploration of the background of the angular momenta dynamics 
leading to specific unison motion for a perfectly entangles qubit pair [37] which is also 
the source of the violation of Bell's inequalities. The main properties of two entangled 
qubits and the results presented in the paper can be summarized as follows: 

(i) Two entangled qubits represented by the quantum state Eq. (CQ) exhibit a 
distribution of angular momenta with projections forming in the x?/-plane an average 
angle (0) = ip. 

(ii) The ensemble average cosine of this angle is to an excellent approximation given 
by (cos(0 — (4>))) ~ sint?. This means that for perfectly entangled states, $ = |, all 
pairs of momenta in the ensemble form an equal angle <fi = if. In general, a higher 
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degree of entanglement is directly related to a pronounced unison motion of angular 
momenta, i.e., with suppressed relative angle fluctuations in agreement with the results 
of a similar treatment in the framework of standard quantum mechanics [12]. 

(hi) In addition to the geometrical picture based on relative angles of the angular 
momenta projections the entanglement of formation can quantitatively be expressed by 
appropriate Bohmian probability distributions and the corresponding Shannon entropy. 

(iv) The virial theorem for a general system of two rotors is confirmed numerically: 
the ensemble averaged quantum potential represents one half of the average kinetic 
energy of the rotors. A general consequence is a known distinction of Bohmian 
pre-measurement values from after-measurement standard quantum mechanics results. 
Specifically, spin-spin correlation tensors differ from customary quantum mechanics 
counterparts by a factor ~ 2/3. 

(v) For fully entangled qubit pairs the Bohmian spin-spin probability distributions 
exhibit power-law tails and are for particular cases given analytically. 

(vi) It is shown that the Bohmian analogue of Bell's inequalities, expressed in 
terms of Bohmian spin-spin correlators, is for fully entangled states identical to the 
standard quantum mechanics counterpart, which demonstrates the non-local nature 
of hidden variables in the present approach. The source of non-locality lies in the 
quantum potential which generates an instant coupling between the angular momenta 
of entangled qubit pairs. In fully entangled states the angular momenta precess in a 
particular manner forming a constant relative azimuthal angle (p which is the origin of a 
specific form of both, the standard quantum mechanical and the Bohmian correlators, 

(a x b x + a y b y ) cosy? + (a y b x - a x b y ) siny? - a z b z , 

leading to the violation of the Bell inequalities. It should be noted that the 
equality of C(a, b) and -B(a, b) guarantees equivalence of Bohmian and standard 
quantum mechanics regarding the (non)violation of all inequalities which can quantum 
mechanically be expressed in terms of the elements of spin-spin correlator tensor 
normalized as in Eq. (T4"T|) . including recently introduced tests for local realistic models 
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